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Abstract 

We calculate the divergences of the generating functional of quenched 
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appropriate definition of the counterterms. We show that the quenched 
chiral logarithms can be accounted for by defining a renormalized -Bo pa- 
rameter which, at lowest order, is proportional to the vacuum expectation 
value of the scalar quark density. Finally, we calculate several quantities 
at one loop to better analyze the modifications induced by quenching in 
the ultraviolet finite part of the one-loop corrections. We point out that 
some of the finite loop corrections may diverge in the chiral limit. 
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1 Introduction 



Most lattice calculations of QCD in its non-perturbative regime and weak in- 
teractions use at present the quenched approximation, i.e. neglect the effect of 
virtual quark loops. Taking them into account considerably increases computing 
times. This means that presumably the quenched approximation will remain with 
us for quite a long time: even with computers much faster than those presently 
available, it will always offer the chance to make a low cost exploratory calculation 
before embarking on a full QCD simulation. 

Simulations of quenched QCD would be much more useful if we had a real 
understanding of the effects of this approximation. Investigations in this direc- 
tion have been made by several authors |1|, 0, j3|. At present we see one main 
approach that has proven to be the most systematic, and also to incorporate most 
of the useful ideas that have been proposed on the subject. This method is called 
quenched Chiral Perturbation Theory (qCHPT), and has been originally pro- 
posed by Bernard and Golterman in Ref. |3j for the purely strong sector (strong 
interactions in the presence of external fields). It has been recently extended to 
the heavy-light meson sector [|J, to vector mesons p| and to the baryon sector 
H . It has been also used in the context of non-leptonic weak interactions @, [?J . 

Let us shortly review the main ideas behind this approach. The difficulty 
to control the quenched approximation comes from the fact that one is modify- 
ing the theory at the non-perturbative level. On the other hand we know that 
at low energy it is possible to define a perturbative scheme to study the strong 
interactions: this scheme is known as Chiral Perturbation Theory (CHPT). In 
this framework the expansion parameter is given by the energy of the weakly in- 
teracting Goldstone bosons of the spontaneously broken chiral symmetry: these 
have a vanishing interaction at zero energy, as symmetry dictates. The chiral 
symmetry imposes also a set of relations between the coefficients of this expan- 
sion in different amplitudes. Those relations do not fully constrain the theory 
that at each order of the expansion has a number of free constants. These con- 
stants incorporate the effect of the non-perturbative QCD dynamics. Under the 
assumption that in the quenched approximation the mechanism of spontaneous 
chiral symmetry breaking is preserved, one may attempt to construct a pertur- 
bative scheme for the quenched case, analogous to the one valid in the full QCD 
case. In this manner one would be able to calculate those effects of quenching 
that modify the perturbative, calculable part of the theory. On the other hand, 
the changes in the unconstrained low energy constants remain unknown, being 
due to the modifications which affect the non-perturbative QCD dynamics. This 
method has the advantage of introducing from the start this clear, useful separa- 
tion between the non-perturbative dynamics of the fundamental theory and the 
perturbative, predictable dynamics of the Goldstone bosons. 

A peculiar aspect of the quenched approximation comes from the U(l) axial 
anomaly of QCD. In the fundamental theory the would-be-Goldstone boson (the 
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rf) does not become massless in the chiral limit, since the axial anomaly generates 
a singlet component (heavy) mass at the level of the effective theory. Thus, in the 
real world the rf is heavy and decoupled from the octet of the pseudo-Goldstone 
bosons. In the quenched approximation this decoupling stops halfway: only one 
of the diagrams that are responsible for the decoupling of the rf survives. At 
the level of the effective theory this has important consequences: the singlet field 
remains light (degenerate with the Goldstone bosons) and has to be treated on 
the same footing as the octet fields. However its two-point function develops 
a double pole and does not admit an interpretation as a propagator. Treating 
the singlet as a dynamical degree of freedom brings in new constants in the 
effective theory. One of them is a new mass scale (the singlet mass m ) that 
is generated by the anomaly, and that does not vanish in the chiral limit. This 
mass appears in the numerator of the double-pole term in the singlet two-point 
function. As different authors have shown [3], this double pole is responsible 
for the presence of a new type of chiral logarithms (we denote them as quenched 
chiral logs) in loop corrections, of the form m^lnM^, as opposed to the standard 
ones M%lnM%. This is one of the main qualitative differences that arises from 
the quenched version of CHPT. 

So far, works in quenched CHPT have concentrated on specific processes, 
analyzing the changes induced in Goldstone boson loops and the size of the effect 
of quenched chiral logarithms. The aim of the present work is to perform a 
complete renormalization of the theory at the one-loop level, on the same line as 
what has been done by Gasser and Leutwyler in the ordinary CHPT case |S|. . 
This requires a calculation of all the ultraviolet divergent pieces of the generating 
functional and a definition of the Lagrangian at order p 4 , the next-to-leading 
order. The advantages of the present analysis are the following: 

1) The calculation of the divergences and renormalization can be done for a 



generic number of flavours N. As we have shown in Ref. |L0[ the N- 
dependence of the divergences can be used to verify the cancellation of 
quark loops in the effective theory. 

2) Like in the standard case, the calculation of the divergences at the generat- 
ing functional level provides a useful check for single amplitude calculations. 
This check is even more welcome in qCHPT where the number of graphs 
to be computed becomes soon very large. 

3) This calculation allows to have full control on the divergences due to singlet 
loops. In particular we will show that quenched chiral logarithms can be 
accounted for via a renormalization of the low-energy constant B (which 
is proportional to the qq condensate). This constant appears in all other 
quantities through the pion mass squared, with the only exception of qq 
matrix elements, that have it as an explicit factor. 
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After having performed the one-loop renormalization, we will devote our at- 
tention to the ultraviolet finite part of the one-loop corrections, by computing 
specific physical quantities at one loop. The relevance of the finite part of the 
loop corrections is in the fact that they may contain terms which diverge in the 
chiral limit like an inverse power of the quark mass. One can realize that this may 
happen by simply looking at the standard chiral power counting [11], and taking 
into account the fact that in quenched CHPT a new vertex appears with chiral 
order zero (the vertex proportional to tuq). Power-like chiral divergences and 
quenched chiral logs are the crucial problem of the quenched version of CHPT: 
the effective theory is defined as an expansion around the chiral limit, and this 
limit is no more well defined in the quenched case. On the other hand these di- 
vergences seem to be unavoidable in the present framework and it looks plausible 
that they are a direct consequence of the sicknesses of quenched QCD. To clarify 
this very important point, a direct evidence of these effects in lattice simulations 
of quenched QCD would be most welcome. 

In our analysis of various observables we will give the complete one-loop 
results. Our aim is not just to make predictions, or to compare with numbers 
produced in lattice simulations. Rather, we would like to show in detail how 
the quenched approximation distorts the matrix elements. For this reason we 
will only work in Minkowski space-time: all the formulae will be given with the 
idea that one should be able to easily see the difference from the corresponding 
ones calculated in standard CHPT. In particular we will stress the presence of 
terms divergent in the chiral limit and of unphysical threshold singularities in 
Minkowski space-time at infinite volume. These type of singularities have been 
already discussed in the literature [|TJ, [TJJ, and have lead to the conclusion that 
quenched CHPT makes sense only in Euclidean space-time. Despite this, we 
still prefer to calculate amplitudes in Minkowski space-time, considering them as 
formal expressions. As we just said this will make the comparison to standard 
CHPT amplitudes easier; on the other hand, the modifications needed to go to 
Euclidean space-time can be easily implemented. 

The plan of the paper is as follows. In Section 2 we outline the main steps from 
CHPT to its quenched version. We give the leading order Lagrangian and define 
our notation, both for CHPT and quenched CHPT. In Section 3 we calculate 
the divergences of qCHPT to one loop using the background field method, while 
Section 4 contains the list of counterterms for a generic number of flavours N 
and for N = 3 and 2. This completes the renormalization of the theory at the 
one-loop level. In Section 5 we analyze a few quantities to one loop in the two 
degenerate flavours case. These are the qq condensate, the pion mass and decay 
constant, the scalar and vector form factors of the pion, and the irn scattering 
amplitude. In Section 6 we state our conclusions. We have also three appendices. 
In Appendix A we give a simple derivation of the divergent term proportional to 
tuq in the quenched generating functional. In Appendix B we give the explicit 
TV-dependence of the divergences in the non-leptonic weak interactions sector, 
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and guess the divergences in the quenched case by simply dropping any N- 
dependence. Finally, in Appendix C we give the explicit expressions for the 
one-loop functions which enter the calculations. 



2 From CHPT to its quenched version 

In this section we introduce the standard notation of Chiral Perturbation Theory, 
that will be also used in its quenched version. We work in Minkowski space-time 
in both cases for ease of comparison. For any further detail in the derivation 
of the CHPT Lagrangian we refer the reader to the original works by Gasser 
and Leutwyler || [|. The construction of the CHPT Lagrangian is based on the 
identification of the symmetry group of the QCD Lagrangian in the chiral limit, 
which, for N flavours is given by U(N) L ® U(N) R , and on the well supported 
assumption that the symmetry of the subgroup SU(N) L <g> SU(N) R is sponta- 
neously broken to SU(N)y. The extension of this construction to the quenched 
case was proposed by Bernard and Golterman || on the basis of an observation 
made by Morel [0]. He observed that, formally, a Lagrangian corresponding to 
quenched QCD can be obtained by adding to the QCD Lagrangian a term which 
is totally analogous to that for quark fields, but which contains ghost spin-1/2 
fields with wrong, i.e. bosonic statistics. The symmetry of the resulting La- 
grangian in the chiral limit is larger than that of QCD and is given by the graded 
group: U(N\N)l ® U(N\N)r, describing transformations between N physical 
flavours and N ghost flavours. It is then assumed that also in the quenched case 
a spontaneous symmetry breaking down to the diagonal subgroup SU(N\N) V 
occurs. Like in standard QCD, the U(1)a symmetry is anomalous. 



Chiral Perturbation Theory describes the dynamics of the octet Goldstone bosons 
fields (pions) of the spontaneously broken chiral symmetry of QCD. It is an 
expansion in powers of the energy of the Goldstone bosons and the light quark 
masses. The lowest order CHPT Lagrangian, i.e. at order p 2 and linear in the 
quark masses, can be written in the following form: 



where (...) stands for the trace over flavour indices, F is the bare pion decay 
constant and the fields are defined as follows 



2.1 Standard CHPT 



C 2 = —(DyUD^U^ + E7+x + X ] U) = —{u^ + X +) , 




(2.1) 



U = u 2 




X 



d„U - irJJ + iUl t 
2B (s + ip) , 
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X+ = u\ ] u + u^u 1 " . (2.2) 

The Lagrangian contains the external sources s,p, f M , a^, = t> M + a M , Z M = 
v^ — a^, which are N x N matrices, with N the number of flavours. The field cf) is 
an Nx N matrix that contains the Goldstone bosons fields: <\> = l/y/2 X^i -1 Xi4> 1 - 
In case, one may add to a singlet component, so that (0) = <p . In the presence 
of a singlet component the Lagrangian in Eq. (|2.1|) is invariant under U(N)l ® 
U(N) R . Since in QCD the U(1)a subgroup is anomalous the breaking pattern 
U(N) L ® U(N) R -> Sf/(iV) L ® SU(N) R ® [7(l) y is realized. The invariance 
under the residual unbroken group allows for the presence of extra functions of 
the singlet component 0o only. A possible choice for this Lagrangian, compatible 
with P, C, T and chiral invariance, is (see also |§ for a different choice) : 

C 2 = F 1 (0o)(^^^ t ) + V r 2 (0o)(^ t X+X t ^)-Vo(0o)+^(0o)^0o^o , (2.3) 
where all the functions are even and real functions of 0o- 



2.2 Quenched CHPT 

The modification needed to construct the quenched version of the CHPT La- 
grangian in Eq. (|2.3|) for a generic number of flavours N consists of the extension 
of the chiral symmetry group SU(N) L £g> SU(N) R ® U(l)y to the graded group 
[SU(N\N) L ® SU(N\N) R ] QU(1) V , which enlarges the spectrum of the theory 
to include ghost states (the © stands for the semidirect product of U(l)y, which 
does not commute with transformations that exchange particles with ghosts). In 
the quenched case there are N physical flavours and N ghost flavours. Under 
the graded extension all the N x N matrices representative of the original U(N) 
group are transformed into graded 2x2 block matrices 

A ^ ( C D ) ' 

whose components are in turn N x N matrices. The matrices A and D {B 
and C) have bosonic (fermionic) character. The trace is then transformed into 
sup ert race: 

tr(A) -> str ^ ^ £ ) = tr( ^ ~~ ' 

The leading order Lagrangian of quenched CHPT can be written in full analogy 
to the standard CHPT caseQ: 

C 2 = V 1 (%)stT(D^U s D^Uj) + V 2 (%)stT(xlU s + Ulxs)-V (%) 

+V 5 (^o)D^ D^ , (2.4) 

3 To distinguish between a quenched CHPT quantity and its standard counterpart we use 
either capital letters (as in <j> — * $) or, when this is not possible, the s subscript (as in U — » U s ). 
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where again Vi($o) are even and real functions of the generalized singlet field $o- 
The graded meson field is defined through the usual exponential representation: 

U s = exp(v / 2?$/F) , 

where F is the bare quenched pion decay constant and $ is now a hermitian non 
traceless 2x2 block matrix 



$ = [l ~ , str($) = $ c 




which contains the new ghost states of the quenched spectrum. All the possible 
quenched meson states carry the quantum numbers of a two particle bound state 
made up with quarks q or ghost-quarks q. On the diagonal sites it contains the 
physical pseudo-Goldstone boson matrix (i.e. the physical pions including the 
singlet component), with the quantum numbers of a qq pair, and the ghost field 
matrix 0, with the quantum numbers of a qq pair. They are both of bosonic 
nature. In the off-diagonal sites are the ghost hybrid fields 9 and 6\ which carry 
the quantum numbers of a mixed qq and qq pair respectively, both of fermionic 
nature. This spectrum of meson states can be found also in the original derivation 
by Morel [l]. He calculated the functional integral over the quark and ghost- 
quark fields (in the leading large-e? expansion and strong gauge coupling limit) 
and obtained exactly the meson spectrum of quenched CHPT, with mesons of 
composite nature, given by bilinears of quarks/ghost-quarks at the same lattice 
site. 

The covariant derivative over the field U s is defined as D^U S = d^U s — 
ir£U s + iU s l%, where r%(l%) is the right (left)-handed external source of the graded 
group. The field Xs = 2B (s s + ip s ) contains the external scalar (s s ) and pseu- 
doscalar (p s ) sources analogously to the ordinary CHPT case. All the exter- 
nal fields r%,l%, s s ,p s are generalizations of the standard external fields, in order 
to make the Lagrangian in Eq. ( |2.4| ) locally invariant under the graded group 
[SUl(N\N) <S> SUr(N\N)} U{l)y. Since we are not interested in studying ma- 
trix elements containing the spurious fields as external legs, we will always use 
the standard external sources only. With this reduction a generic graded source 
reads as follows: 

3s = [ j ' 3 = v ' v ^ a ^ ' 

For the scalar external source we must recall that it is defined to contain the 
quark mass matrix Ai which has to be the same both for the quarks and the 
ghosts: 

( M + 5s \ 
Ss ~ { M ) ■ 

In what follows the quark mass matrix will be taken proportional to the unit 
matrix: M. = m q l. All the Goldstone bosons and their ghost counterparts will 
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have the same mass: M 2 = 2B m q . We have adopted the usual CHPT notation 
and call M 2 the lowest order term in the expansion of the mass of the pions in 
powers of quark masses: 

M 2 = 2B m q + 0{m 2 q ) = M 2 + 0{m 2 q ) . (2.5) 

Finally, we expand the functions V^($ ) m powers of $ : 

Vi,2(*o) = ^ + \^l + 0{^) , 

f)' 

V 5 (%) = We +0(*&) , (2.6) 
and we shall always work with number of colours N c = 3. 



3 One— loop divergences 

To calculate the ultraviolet divergent part of the quenched generating functional 
to one loop we use the background field method, i.e. expand the action around 
the classical solution, which is determined by the external sources through the 
classical equations of motion. We write the field U s as: 

U s = u s e tr - u s , 

where U s = u 2 is the classical solution to the equations of motion. In the absence 
of spurious external sources it reduces to 

/ u \ 

Us = { o i J • 

We decompose the fluctuation S similarly to the field $ and write: 

E =[l £ ) > str(S) = vOvfo - 60 , 

(note that with this normalization the £o an d Co have a proper kinetic term). The 
matrix fields £ and ( are decomposed as follows 

N 2 -l N 2 -l 

e= E u 4 , c= E **c, (3.i) 

a=0 a=0 

and the fields £ and analogously, where A a = \ a /\/2, a = 1, . . . , N 2 — 1, and 
A = 1/VN. Given their special character, it is useful to separate the singlet 
components of the £ and £ fields from the rest, and combine them into one vector: 

Mi)- 
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The remaining fields are put into components of the following vectors: 



? = te, e,...,^ 2 - 1 ), c f = (c t0 , c n ,..., c^ 2 - 1 ). 

With this notation the action can be written as 

Sl$] = S[$] ~^Jdx {X^D X X + + £oB T Z + fB£ + 2C f A;C 

- i T (D + M 2 )i} + 0(E 3 ) . (3.2) 

The explicit expressions for the various differential operators Dx,(£ and the ma- 
trix B will be given below. The matrix B induces a mixing between the singlet 
and non singlet component of the physical meson field. Notice also that the fields 
£ are completely decoupled from the rest: the integration over these degrees of 
freedom produces only an irrelevant constant. 

Before deriving the various contributions to the generating functional at one 
loop we shift the field £ in order to remove the mixing with the singlet component 
£o. By performing the translation 

£ = £' - D^B£ , 

one gets 

£ T D£ + £ 5 T £ + fB£ = ^D^' - ZoBTD^BZq . (3.3) 

In this manner the action up to the quadratic fluctuations becomes a sum of 
quadratic differential forms diagonal in the fields X , £, £, £, £t. The price 
to pay is that now the differential operator acting on the singlet field Xq has a 
nonlocal term. Denoting as Dx the new non local operator acting on Xo after 
the shift, the quenched generating functional to one loop can be formally written 
as follows 

iz^HPT det£> c 

e noop = Al 1 r • 3.4 

(det J D f )2(det J D x )2 

As we will see below, the non locality of Dx will hardly make the calculation of 
the divergent part more complicated. 

3.1 Integral over the £ fields 

The differential operator D^ b is defined as follows^ 



d,C = d^ a + t% , (3.5) 



4 We remind the reader that the indices a, b run from 1 to TV 2 — 1. The singlet components 
are treated separately. 
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where 



-(r,[A a ,A 6 ]) 



«■ =-\(K,H^^]) + \(& a A b }x + ) , (3-6) 



and = l/2([u*, d^u] — iu^r^u — iul^v)) is the vector current connection of the 
covariant derivative over the dynamical fields. 

The ultraviolet divergent part of the integral over the £ fields can be derived in 
closed form by regularizing the determinant in d dimensions and using standard 
heat-kernel techniques. The result reads: 



In det Dp 



+ 



(4vr) 2 (d-4) 
■(( 

N 1 \ 



+ g-<(«X) 2 ) + j(«/x+) + \(uX)(x+) 



8 2NJ 



8 47V 2 
+ . . 



(x + y 



(3.7) 



where the ellipsis stands for contributions which are finite in four dimensions. 
This result is the standard CHPT result derived in ||, where now we also keep 
terms proportional to (u^) that are nonzero only in the presence of the singlet 
component. 



3.2 Integral over the £ fields 

The differential operator D^ b is defined like in Eq. (|3.5|) , but with barred quan- 
tities, given byQ 



Ttt6 



(T,X a X b ) , a c 



|((«/ + X + + ^M)\ a \ h ) , (3i 



where we recall that *M is the quark mass matrix. 

The ultraviolet divergent part of the functional integral over the ( fields can 
also be given in closed form using standard heat-kernel techniques. The result 
reads: 



i In det D 



-1 



C 



(4vr) 2 (d-4) 



dx 



jiT^rn + ^((n/ + X+ + 45o-M) 2 ) 



+ . 



(3.9) 

As we remarked in Ref. [1C], the integral over the ( fields completely removes 
the terms linear in N in the divergences of standard CHPT to one loop. This 

5 Here the singlet component is included, and the indices a, b run from to N 2 — 1. 
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dependence is not fully explicit in Eq. ( |3.7| ), 
trace of x+> when we expand this around s ■ 



since a factor N is contained in the 
= A4 and for A4 diagonal: 



{ X+ ) = 2NM 2 + 0(<P 2 ) . 

This result shows that the qCHPT scheme is coherent: the terms linear in N can 
only be generated by quark loops, and these are supposed to be absent in the 
quenched approximation. 



3.3 Integral over the Xq fields 

After the shift of the £ field the operator acting on X can be written as: 



where 



D 



x 



u x 



A x 



D x --(l + r 3 )B T D^B 



(3.10) 



D° x + A 



x 



N 



r 3 (n + M ) + — (1 - T X ){aU + m, 



4N 



(i + r 3 )(x+) - N(i - n ) {v x {uX) + Mx + )) + o($o) 



(3.11) 



and x + = X+~ 2M 2 1, so that (%+) = (x+) ~ 2NM 2 . The expression of D° x , 
the "free" part of the differential operator, clearly shows that the theory has a 
problem here: it is not possible to diagonalize that operator, and we do not have 
two freely propagating normal fields (£o>£o)- O n the other hand this problem is 
welcome in this context, since it is thought to be the manifestation of the absence 
of quark loops in the singlet field propagator, at the level of the effective theory. 
In the language of Feynman diagrams this problem shows up as a double pole in 
the propagator of the singlet field, whose consequences on observables have been 
studied by several authors p], ||, ^2| . We adopt the usual point of view on this 
problem, i.e. assume that it has to be there, and proceed with the calculation of 
the divergent part of the generating functional. 

In this case we cannot apply straightforwardly the heat-kernel techniques, 
because the differential operator does not reduce to a diagonal Klein-Gordon 
operator when the external fields are put to zero. Therefore we just expand the 
logarithm of the differential operator, and isolate the ultraviolet divergent terms: 



Trln (D x /D° x 



Tr 

iiv 

2 



D X -\D X -D 



x, 



D X -\D X -D X )D X -\D X -D X ) +... .(3.12) 
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One can easily see that the ellipsis in ( |3.12| ) contains ultraviolet finite terms only. 
We postpone a more detailed discussion of the infrared behaviour of Eq. (|3.12| ) 
to the end of this section. The inverse of the "free" operator D x is: 



D 



o -i 

X 



Gri 



r 3 -(l+r 1 )y(an + m2)G 



where 
and 



(n + M 2 ) x G (x-y) = 6(x-y) 



D x -D 



x 



A x --(l + r 3 )B T D7 1 B . 



(3.13) 

(3.14) 
(3.15) 



As we anticipated above, the overall effect of the shift made to remove the mixing 
between singlet and nonsinglet fields is easily accounted for. Expanding around 
the free part of D7 1 in the non local term one gets 



B T D7 l B 



AN 



+ 0{Gl) 



(3.16) 



The term proportional to 0(Gq) can only yield ultraviolet finite contributions 
to Eq. ( |3.12| ), while the Go term yields ultraviolet divergent contributions only 
to the first term of the expansion in Eq. ( |3.12| ). This shows that also in the 
singlet sector the UV divergent part is local and can be given in closed form. 
The calculation of the ultraviolet divergent part of In det Dx is now easy: we 
simply have to insert back Eqs. (|3.13|J3.15|) into Eq. ( |3.12|) and keep only the UV 
divergent parts. Having worked out the traces (over the r matrices) we obtain: 



Trln (d x /D 



(4vr) 2 (d-4) VAN 



(X+) gjy2^+) 



+ - 



6 
1 



- 7,(x+) M«X> + v 2(x+)) 



+ 



(3.17) 



where the ellipsis contains UV finite terms only. The terms proportional to inverse 
powers of N exactly cancel those contained in Eq. fl3.7| ) giving a result that is 
totally iV-independent. The terms proportional to m\ and powers of a are the 
effect of the double pole in the singlet propagator, and are also iV-independent. 
Note that no mixed terms of the type (rriQ,a) x (fi,f2) can be produced in the 
divergent part. 

The term proportional to rrv^ is a term already present in the 0(p 2 ) La- 
grangian. To remove that divergence one has to add to the lowest order param- 
eter B in the £ 2 Lagrangian a <i-dependent part proportional to Wq that has a 
pole at d = 4: 

1 2ml 



B t) 



Bn 



d-A 



1 + 



16tt 2 d - 4 3F 2 



(3.18) 
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This feature is completely new with respect to standard CHPT (in dimensional 
regularization), and stems from the fact that in the quenched theory we have 
a new mass scale that does not vanish in the chiral limit. After the divergence 
has been removed, we are left with a term of the form m 2 ) l i iiM 2 (x+)- This term 
contains all the one-loop quenched chiral logs that have been discussed at length 
in the literature. Our calculation shows that they can be fully accounted for by 
defining a renormalized constant Bq 

B o^B = B ^l-^ r2 \n— + bM ] j . (3.19) 

Notice that since Bq is independent from the quark masses, Bq becomes divergent 
in the chiral limit. To find evidence for these quenched chiral logs one should try 
to extract from lattice data this quantity Bq. As we will see the quark condensate 
and the scalar form factor are two excellent candidates for this, since they are the 
simplest quantities which are explicitly proportional to Bq. Other quantities will 
tipically depend on Bq through the renormalized pion mass. This at one loop is 
given by: 

M 2 = 2BQm q + 0(m 2 q ) , (3.20) 

and is not divergent in the chiral limit. These other quantities are therefore much 
less suitable to identify the presence of quenched chiral logs. 

Of course what we have just said is valid in the specific sector we are studying 
here. To extend it to other sectors of the effective theory (like the non-leptonic 
weak interactions) requires further study. However we have a rather simple argu- 
ment that shows that what has happened here will happen also in other sectors: 
the quenched chiral logs to one loop contribute to the redefinition of one of the 
constants appearing in the lowest order Lagrangian. In order not to interrupt 
the discussion here we relegate the argument to Appendix [A]. 

3.4 Complete result 

In this section we put together all the various pieces and give the complete result 
for the ultraviolet divergent part of the generating functional of qCHPT to one 
loop. The explicit expression for Eq. (|3.4| ) is: 



7 qCHPT \ 1 J 

~ (4vr) 2 (rf-4) ' ' 



1 loop 



l -(u,u v )(u^) + ^(uX) 2 



+ - (1 - Av x ) + -a (! - 8 ^) (x+Y 



I 1 

16 

2 2 

+^(x+) + ^(x+) 2 -^(x+) 



7W(«" {u v u" + x+)) 



+ ... . (3.21) 
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The most striking feature of Eq. ( |3.21 ) is the complete flavour independence of 



the result. If we analyze in detail the modifications that the quenched approx- 
imation has produced to the divergent structure of the effective theory at the 
one-loop level, we come to the following list: 

1. all the terms proportional to N have been dropped; 

2. all the terms proportional to 1/N and 1/N 2 have been dropped; 

3. new divergences proportional to the parameters present in the anomalous 
singlet sector have been produced. 

All these new parameters are dimensionless, with the only exception of itiq. The 
dimensionless parameters (a and v\^) generate divergences that have the struc- 
ture of a chiral invariant term (since they do not break the chiral symmetry) of 
order p A , for obvious dimensional reasons. For the same reasons generates 
divergences with the structure of a chiral invariant of order p 2 . As it is shown in 
Appendix |A| one can very easily understand why it is only the mass term (x+) 
that is generated. 

As it turns out, the modifications listed in points 1. to 3. above find a very 
simple explanation: dropping the terms proportional to iV corresponds to drop- 
ping virtual quark loops. Dropping the terms proportional to 1/N and 1/N 2 , is 
a consequence of having a singlet degenerate in mass with the nonsinglet pseu- 
doscalars. The new parameters in the singlet sector are required by the U(1)a 
anomaly, and the diseases in that sector are generated by the quenched approx- 
imation, as it is well known. These simple conclusions suggest that one could 
have guessed all these modifications without doing any calculation. In fact, we 
provide an example of how one could try such a guess in Appendix |B|, where 
we apply the same criteria to the generating functional of the non-leptonic weak 
interaction sector for the octet on-shell case (the complete analysis will be given 
elsewhere by going through the three steps we have enumerated above. 



3.5 Chiral and threshold divergences 

Quenched chiral logs are not the only problem generated by the presence of the 
double pole in the quenched version of the singlet propagator. As we will see 
in detail in Sec. 5 through several examples, this double pole generates also 
other kind of divergences inside contributions which are ultraviolet finite. These 
divergences are of two types: powerlike chiral divergences, i.e. inverse powers of 
M 2 , and unphysical threshold divergences. We find it instructive, before closing 
this section to identify which are the terms in the generating functional which 
are responsible for them. 

Some of the terms (and in fact an infinite series of them) that we have ne- 
glected in Eq. (|3.17|) because they are ultraviolet finite, contain these kind of 
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singularities. They can be given in closed form only if one stops at a given order 
in the expansion in powers of the field $. Since in the following sections we are 
not going to analyze anything beyond the four-point function, we can stop at 
order $ 4 , and identify explicitly the troublesome terms. They all come from the 
insertion of the double pole term of ( |3.13| ) in the expansion ( |3.12| ), and give the 
following contribution to the generating functional: 

SZfZ T = K ~ 4 QM2) / dxdy h{x - y)(x + (y)x + (x)) 

- a J* M - / dxdy I x (x - y)(x + (y)){x + (x)) 



72 

(mg - aM 2 ) 2 
144 



dxdy I 2 {x - y)( X+ {y))(x + {x)) + 0{§ ( ) .(3.22) 



The functions h(z), l2(z) are defined in appendix |CJ. At infinite volume and in 
Minkowski space-time their Fourier transforms I\(q 2 ) 1 l2(q 2 ) develop an imagi- 
nary part when q 2 > AM 2 which diverges at q 2 = 4M 2 (see appendix 0). More- 
over their values at q 2 = are inversely proportional to M 2 (again see appendix 
0): this is the origin of powerlike chiral divergences that we will find in several 
observables in Sect. 5. The threshold singularities in particular make the theory 
meaningless in Minkowski space-time at infinite volume. In finite volume and in 
Euclidean space-time the same one loop functions I\^{q) have been evaluated at 



q = AM 2 in Ref. [12], and it was found that these functions give rise to enhanced 



finite volume corrections which are forbidden in a healthy Hamiltonian theory. 



As it was pointed out in Ref. [12] this shows that qCHPT can only make sense 



in Euclidean space-time and in finite volume. 



4 Lagrangian at order p 

To complete the renormalization of the quenched theory at order p 4 one needs to 
add the most general chiral invariant Lagrangian at this order. As in the standard 
CHPT case, some of the couplings appearing in the order p A Lagrangian have an 
UV divergent part in such a way that all the one-loop divergences are removed. 
The most general chiral invariant Lagrangian at order p A in standard CHPT has 
been given by Gasser and Leutwyler ||. The extension to the graded symmetry 
version is not needed here, since we are not going beyond order p 4 , and are not 
interested in having the spurious degrees of freedom as external particles^]: we 
can use the standard CHPT Lagrangian right away. 

There is however a slight modification that we have to introduce. As we 
noted before, the trace of x+ starts with a constant term proportional to N in 

6 Moreover we will not consider singlet fields as external particles. They require at least two 
more counterterms as shown by Eq. ( |3.21 ). 
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the degenerate mass case we are considering here. In the quenched version a 
linear dependence upon N is forbidden, and therefore we must always substitute 
(x+) —> (x+)- Apart from this modification, we have followed existing notations 
for the choice of the 0(p ) Lagrangian, both in the SU(3) and £77(2) case. The 
SU(3) choice is the standard Gasser and Leutwyler Lagrangian H], while for 



SU(2) we choose to use the Gasser-Sainio-Svarc Lagrangian fT4|j . 

An important point concerns the value of the counterterms: we observe that in 
the quenched case the counterterms do not depend on the number of flavours. Not 
only the divergent part, as we have explicitly shown in the previous section, but 
also the numerical value of the finite part of the counterterms does not change 
for different values of N. Therefore it is useful to identify, and give names to 
them in the general N case. For the more interesting cases of N = 3 and N = 2, 
because of trace relations, one will be able to access only certain combinations of 
them, as we will specify below. For general N the lagrangian at order p 4 is given 
by: 

10 

£ 4 = £A?P, , (4.1) 

where the eleven operators Pi are listed in Table [l] (we remind the reader that in 
the quenched case it is necessary to change (x+) — » (x+))- These eleven chiral 
invariant operators contain, besides those defined in Eq. ( |2.2| ), the following new 
building blocks: 

X- = u ] xu ] - ux ] u . (4.2) 
To derive the results shown in Table [l] the following relation is useful: 

f +flu = 2i r Mi , - -[up, Uu] , (4.3) 

and the identification {%%) = 1/2 (x+ + X-) and (/|) = 1/2 (/£ - f_) can be 
done up to contact terms which contain external sources only. 

4.1 SU(3) Lagrangian at order p 4 

For N = 3 the Lagrangian at order p A reads as follows: 

10 



C[ N = z) = Y.m , (4-4) 
1=1 



where the operators Pi are defined in Table |^. The Pq operator is linearly depen- 
dent from the others through the following trace relation: 

Po = \Pi + P2- 2P 3 , (4.5) 
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Table 1: List of terms of order p 4 for tV generic, N = 3 and N = 2. In the 
second and third columns we give the coefficient of the divergence coming from 
the one loop in the standard and quenched CHPT case. As we have indicated in 
the table, the invariants containing (x+) have to be changed with (x+) —> (x+) 
in the quenched case. 



16 



which implies: 



(4.6) 



In order to reabsorb the divergences at one loop we define the Lf in the following 
manner: 



,.d-4 r 

A 



16tt 2 



_l__i (ln4 . + r'(i) + i) 



(4.7) 



with ji the renormalization scale, A contains the divergence at d — 4 and the 
coefficients V- are given by 



r? = 



16 ' 
1 

16 



rg = ^ ( 1 - 8u 2 + W 



2- g , 
2 
9 



n = i(i-4 Vl ) , 

r - a 

Fs " "12 ' 



all the other r? are zero. 



(4.8) 



4.2 SU(2) Lagrangian at order p 4 

For N = 2 the Lagrangian at order p 4 reads as follows: 



where 



7 

i=l 



£(iV=2) 



Q 7 



1 

4 



16 



(X- 



(4.9) 



(4.10) 
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To reduce the number of chiral invariants needed we have used the following 
relations: 



Pn 



Pa 



5* 



\ (fl> + Pi) 



(4.11) 



which imply the following relations between the N 



l\ + -a 

1 2 3 



4 l 2 



LI 



3 and N 
1 „ 



16 



T q -4- " T 



2 counterterms: 
1 



~-l q 



Ll + 



±11 
16 L 7 



L% 



2-^8 



(4.12) 



Note that the trace relations have been written down using the invariant (x+), 
and must be reexpressed in terms of i n the quenched case. This generates a 
correction to the constants appearing in the £2 Lagrangian, see below. In order 
to reabsorb the divergences at one loop we define the l\ in the following manner: 



(4.13) 



with: 



it 



1 

4 



7l 



7l 



l-4ui 



7? 



7l 



-a 



1 - 



i% a 

3 



-a 



(4.14) 



all the other 7? are zero. We find useful for the analysis of the phenomenology 
to introduce the scale independent constants 7J, defined as follows: 



32tt 2 M 2 
—5-/? /i - In — 

7i M 



(4.15) 



As we mentioned above, a complete renormalization at the one-loop level re- 
quires, in the quenched case, the renormalization of the order p 2 constant B 
due to divergences proportional to m 2 ,. In the present case, (x+) has been elim- 
inated with the use of the Cayley-Hamilton relations ( |4.11| ) in favour of (x+) 2 
and M 2 (x+). The divergence proportional to the latter can also be reabsorbed 
in the renormalization of the B parameter. Since P5 = 4Q 4 + 2M 2 (u fl u fl ) the 
constant F 2 receives a finite correction proportional to L\. For later convenience, 
we define here the renormalized constants at order p 2 in the two-flavour case, in 
such a way that they include also finite corrections: 

F 2 

■(u fl u fi + x+) , 



Co 



r>JV=2 



B, 



N=2 





1 + 4L?^ 



5 F 2 J 



m: 



2 - 2«M 2 ) 



48tt 2 F 2 



In 



M 2 



8L1 + 



a 



48tt j 



M 2 

+ h M 



(4.16) 
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where, with an obvious notation x+ stands for the analogous of x+ which contains 
B instead of B . 



5 Analysis of various observables in quenched 
CHPT 

In this section we make a complete one-loop analysis of several observables in 
quenched CHPT. The main reason for this is to study the problems generated by 
quenching in the finite part of the one-loop corrections, which we have not con- 
sidered in the generating functional. As we will see, some of the finite corrections 
diverge in the chiral limit. The origin of these divergences can be traced back 
to the presence of the double pole in the singlet two-point function. The double 
pole carries in the numerator a new mass scale mo that does not vanish in the 
chiral limit, and hence modifies the chiral power counting valid in CHPT. The 
standard power counting goes as follows: the chiral order of a generic diagram is 
given by the simple formula 

D g = 4L-2I + J2dN d , (5.1) 

d 

where D g is the chiral dimension of a graph g that has L loops, / internal lines, 
and Nd vertices of chiral dimension d. The topological relation 

L = I-V + 1 , (5.2) 

where V = J2d Nd is the total number of vertices, can be used to obtain 

D g = 2L + Y / (d-2)N d + 2 . (5.3) 

d 

Since in standard CHPT the lowest chiral dimension of a vertex is two, the chiral 
dimension of a graph is always bigger than two, and increases with the number of 
loops and vertices with chiral dimension bigger than two. In quenched CHPT the 
situation changes, and we have to allow for the presence of vertices with chiral 
order zero, i.e. the insertions on the singlet propagators proportional to ??1q (that 
is a constant in the chiral limit). In this case Eqs. ( |5.ip through ( |5.3| ) are still 
valid, but due to the presence of terms with d — 0, D g may now be smaller 
than two, and even negative. Naively one could conclude that D g could even 
be unbounded from below. However, one has to take into account the fact that 
virtual quark loops are forbidden: this puts a series of constraints on the type 
of graphs with 7Bq insertions that are allowed. For example: it is not possible 
for two tUq vertices to lie on the same line one after the other, or, no standard 
vertices can have all the outgoing lines that end up on an vertexj] These 

7 There is one exception to this, given by vertices with physical external sources. In this case 
disconnected quark loops are allowed, since they are not generated by the QCD determinant 



(see Section 5.2) 
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constraints are such that D g comes out to be bounded from below, although it 
may be negative. The value of the lower bound depends on the observable - we 
will see explicit examples below. 

In what follows we are going to analyze: the quark condensate, the pion mass 
and decay constant, the vector and scalar form factors of the pion and the titt 
scattering amplitude. Although these quantities (with the exception of the form 



factors) were already analyzed at the one-loop level in previous works P, g, 0, 
we find it useful to present them here again, in view of the renormalization that 
we have performed at the level of the generating functional, and also of our 
definition of the lagrangian at order p A . We make the analysis in the case of two 
light flavours with degenerate masses. 

5.1 Quark condensate, pion mass and decay constant 

As anticipated in the previous section the renormalized scalar quark condensate 
plays a crucial role among the quenched observables in the strong sector, since it 
contains an explicit dependence upon the quenched chiral logarithms through the 
Bq parameter Q4.16D (everywhere in this section we shall use the B parameter 
as defined in Eq. (|4.16|) , dropping the N = 2 superscript). We shall see later in 
the case of the scalar form factor that all the qq matrix elements share the same 
feature. The renormalized scalar density to one loop in the two-flavour case is 
given by 

(qq) q = -F*Bo[l + 0(M 2 )] , (5.4) 

where we have not written down explicitly the standard chiral corrections of 
order M 2 . The problem with these corrections is that they contain contributions 
coming from counterterms of order p 4 that contain only external fields (we have 
not written them down in the previous section). These counterterms cannot 
be determined on a phenomenological basis: their presence in the expression 
of the quark condensate reflects the fact that away from the chiral limit, this 
quantity cannot be defined unambiguously. We refer the reader to Ref. |§ for a 
detailed discussion of this point. On the other hand, in the chiral limit, where 
this ambiguity disappears, the quark condensate diverges due to the quenched 
chiral logarithms inside B . 

The pion decay constant to one loop is renormalized only by a finite amount 
in the quenched two-flavour case: F v = F, see Eq. ( 4.16[ ). Notice that in the 



quenched three-flavour case there is no need to define an F, but on the other 
hand L\ directly contributes to F n in such a way that for N = 3 and N = 2 
(as also for any other N) one has the same pion decay constant, as expected^. 
The diagrams which renormalize the pion mass to one loop are shown in Fig. 
[I]. The meson tadpole and its ghost counterpart cancel each other, so that the 

8 We thank Martin Liischer for pointing out an inconsistency on this point in the first version 
of the manuscript. 
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Figure 1: One-loop diagrams in quenched CHPT that contribute to the squared 
pion mass M%. They are the meson tadpole, its ghost counterpart and the tadpole 
with one singlet vertex (x) insertion. 

renormalization of the quenched pion mass at one loop is provided by the tadpole 
with one singlet vertex insertion and its counterterm: 

Ml = 2B m q , (5.5) 

where m q is the light quark mass. As one can see, all the one-loop correc- 
tions, including the quenched chiral logarithm have been reabsorbed in B Q . Since 
Bo m q ~ m q logm g when approaching the chiral limit, the renormalized pion mass 
tends to zero like m q \ogm q . No divergence is produced by quenching in the be- 
haviour of the renormalized pion mass in the chiral limit, although the way it 
approaches zero is different from that of standard CHPT. Once M 2 is fixed to its 
physical value no residual quenched chiral logarithms are left in the strong sector 
(with the mentioned exception of qq matrix elements). In Appendix || it is shown 
that the same situation occurs in the weak AI = 1/2 sector, where additional 
quenched chiral logarithms can be reabsorbed in the renormalization of the weak 
mass term. 

5.2 Scalar form factor 

The scalar form factor of the pion is defined by the matrix element of the qq 
density between two pion states 

(7tV)Iw|t*(p)> = S*Fs(t) , (5.6) 

where t = (p — p') 2 . In quenched CHPT the complete list of one-loop diagrams 
which give contribution to Fg(t) are shown in Fig. 0. 

An explicit calculation shows that the fermionic ghost loops do not fully cancel 
the corresponding meson loop diagrams. The reason for this mismatch is best 
understood within the quark-flow diagram picture. Here, the physical scalar 
source only couples to the quark lines and not to the ghost lines. The possible one- 
loop diagrams are the ones listed in Fig. Diagram (b), where the scalar source is 
coupled to the internal disconnected closed quark line, has no correspondent ghost 
loop diagram: this is correct because the loop is not produced by the fermionic 
determinant, and must therefore be present also in the quenched approximation. 
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V 2 




Figure 2: One-loop diagrams in quenched CHPT which contribute to Fg(t) (the 
box stands for the scalar source insertion). . They are the "standard" meson loop 
diagrams (first line) to which one has to add the corresponding fermionic ghost 
loop diagrams, the singlet insertion diagrams (second line) and the diagrams with 
vx, V2 vertex insertions (third line). 



The complete renormalized quenched scalar form factor can be written as 
follows 



1 



Fl 



+l^Ut){rnl-aMlf\ + Oif 



(5.7) 



where the coefficients 7? have been defined in Eq. ( |4.14j ) and Fg(0) is given by 



F q s (0) 



2B <1 



+ 



Ml 



lQn 2 F2 



48vr 2 F2 
7l {%-!) -l\(J\-l) 



2 \ l {ml-aMlY 

3 7 9 48tt 2 F2M2 



(5.8) 



Note that we are working in the degenerate mass case, so that no isospin breaking 
effect has been taken into account. In standard CHPT there is no isospin breaking 
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(a) (b) (c) 



Figure 3: One-loop diagrams which contribute to F$(t) in the quark-flow diagram 
picture. Diagram (b) can be present in the quenched approximation, while the 
others disappear. 

correction to the scalar form factor at this order of the expansion. In passing, 
we state that also in the quenched case there is no isospin breaking contribution 
linear in m u — to the pion scalar form factor, as it happens in CHPT, while 
an isospin breaking correction of order [m u — m^) 2 is produced via the (0o , 
mixing for neutral pions by the chiral invariant operator P 7 in Table [If]. 

The functions J(t), I\{t) and ^{t) are finite and they are defined in Appendix 
[CJ The two functions Ii(t) and hit) are peculiar of quenched CHPT. They will 
also appear in the irn scattering amplitude, where we shall analyze in some details 
the various sicknesses of which they suffer. Here we used their low momentum 
expansion to define the scalar form factor at t — 0. 

The scalar form factor is a good example to analyze the modifications pro- 
duced by the quenched approximation to an observable at the one loop level. 
First, the pion loops have been only partially cancelled, therefore the ordinary 
chiral logarithms and the one-loop function J(t) do appear in the same way as 
in standard CHPT, but with different coefficients (these coefficients may even 
vanish in particular cases, like M w and F n ). Second, quenched chiral logarithms 
appear at one loop, but they can be reabsorbed in the renormalization of the B 
parameter, as we have demonstrated in the previous section. Besides quenched 
chiral logs, the remaining finite loop corrections arising from the anomalous sin- 
glet sector and proportional to are even more problematic, since they have 
negative chiral dimension, as anticipated in the general discussion above. It is 
a simple exercise to calculate the chiral dimension of the one loop diagram with 
two TOq insertions on the two internal singlet lines (this is the central graph in 
the second line of Fig. 0): with respect to the tree level graph this has chiral 
dimension —2. These corrections diverge in the chiral limit like an inverse pion 
mass squared, see Eq. (|5.8| ). In fact, there is an infinite series of graphs that has 
the same chiral dimension: these graphs are obtained from this one by adding 
any even number of singlet lines (each one with one insertion) between the 
two vertices. Also the insertion of tadpoles and sunset diagrams with the maxi- 
mum allowed number of Wq insertions does not change the chiral dimension of the 

9 Note that also the neutral pion mass M 2 gets next-to-leading corrections of order 
O {(m u — to^) 2 ) from Pj in the quenched case. 
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starting diagram. As far as we could see this series of graphs is also the one with 
the lowest chiral dimension for the scalar form factor. This example shows that 
despite the general formula ( |5.3j ) with d = vertices allowed, in the quenched 
case the chiral dimension of amplitudes is bounded from below. 

It is also interesting to look at the slope of the scalar form factor at low 
momenta in the quenched case. This defines the scalar radius as follows 



F§(0) 



(5.9) 



The scalar radius in the quenched approximation at one loop is given by: 
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(5.10) 



In standard CHPT the scalar radius diverges in the chiral limit because of the 
presence of t dependent chiral logarithms. It behaves like: 



lnMi + ... . 



(5.11) 



In the chiral limit the one-loop contribution to the quenched scalar radius di- 
verges not just logarithmically as in the standard case, but like an inverse power 
of the pion mass: 
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+ ... . (5.12) 



The origin of this power-like divergence in the chiral limit is the same as that of 
the form factor at t — 0. Here it is more severe simply because the definition of 
the radius implies a derivative with respect to t. 

It is interesting to note that in quenched CHPT the Feynman-Hellman the- 
orem || |15| does not hold: 

dM 2 

W» * s£ • (5.13) 

as one can easily verify by comparing Eq. (|5.5|) and Eq. ( |5.8|) . The origin of 
the violation of this theorem is in the presence of diagram (b) of Fig. [3] in the 
quenched scalar form factor. This graph cannot be obtained taking a derivative 
with respect to m q of M 2 , since the quark loop is not present in M 2 and cannot 
be resurrected by a derivative. 
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Figure 4: One-loop diagrams in quenched CHPT that contribute to F v (t) (the 
box stands for the vector source insertion). They are the "standard" meson loop 
diagrams (first line) and the corresponding fermionic ghost loop diagrams (second 
line). No singlet component can run in the loop. 



5.3 Vector form factor 

The vector form factor of the pion is defined in terms of the matrix element of 
the vector current = qj^q between two pion states: 

{Ap')KW{p)) q = ie ikl (p,+ P ',)F v (t) , (5.14) 

where t = (p—p') 2 - The divergent contributions to F v (t) can be derived from the 
expression (p. 21 ) of the quenched generating functional in the usual way. It is an 



easy exercise to show that these contributions are zero. In fact, the only chiral 
invariant which can give corrections at order p 4 is the operator number 9 in the 
list of Table [I], which has no divergent term in the quenched limit. In a Feynman 
diagram approach the graphs which contribute to one loop are shown in Fig. f|. 
The complete calculation gives zero, because of the systematic cancellation of 
each pion loop with the corresponding ghost loop. In addition, since no singlet 
component can run in the loop, there is no extra contribution coming from the 
anomalous singlet sector. The quenched vector form factor for N = 2 can be 
written as follows 

F v {t) = l- l f 2 t + 0{p') , (5.15) 



where the finite counterterm l\ is defined in Table p] and Eq. (|4.10| ). Again 



no isospin breaking effects have been taken into account. In standard CHPT 



the Ademollo-Gatto theorem |16] guarantees that they are absent at this order. 
In quenched CHPT the theorem is also valid. Note that the counterterm P 7 
cannot contribute at all to the vector current matrix element, while the new 
chiral invariant term (u fl ){u tl x+) induced by the dynamical singlet component 
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Figure 5: Two-loop diagrams in quenched CHPT with the m 2 , singlet vertex 
(x) insertions that give divergent contributions to the e.m. charge radius in the 
chiral limit. They are tadpoles, which generate quenched chiral logs and the fish 
diagrams that also generate power-like divergences. For each diagram the chiral 
dimension respect to the tree level is given. 

gives O ((m u — m^) 2 ) corrections to the decay amplitude tc + — > 7r°ez/ via the 
(0o , h) mixing. 

Since the vector form factor does not receive contributions from singlets run- 
ning inside the loop at the one-loop level, it does not show any divergence in the 
chiral limit. The situation however changes at two loops, where we have among 
others the graphs shown in Fig. [5| The most dangerous graph is the fish diagram 
with two wig insertions (the last of Fig. |^) which has chiral dimension zero respect 
to the tree level. Again this is only the first example of a full series of graphs 
which have the same chiral dimension: they are obtained from the starting one 
by inserting any even number of singlet lines between the same two vertices as 
those of the two-loop fish diagram, or tadpoles and sunset diagrams all with 
the maximum allowed number of m 2 , insertions. In this case there are no graphs 
which are more singular than those in the chiral limit. 

The low energy representation of F v (t) also determines the electromagnetic 
charge radius of the pion in the quenched approximation 

F v (t) = l + Ur 2 ) v + 0(t 2 ) . (5.16) 
6 

In standard CHPT the presence of t dependent chiral logarithms makes the elec- 
tromagnetic charge radius diverge in the chiral limit || 

(rV = - I ^logM 2 + ... . (5.17) 

The divergence of the electromagnetic charge radius in full QCD can be under- 
stood in a physically intuitive way. The charge distribution is cut off by the 
Yukawa potential ~ e ~ M ^ r at large distances. In the chiral limit M n goes to zero 
and the Yukawa potential is no more effective, the charge distribution falls off 
like a power of the distance and the charge radius becomes infinite. The charge 
distribution of the pion cloud surrounding any particle gets modified by quench- 
ing. As a consequence, the behaviour of the charge radius in the chiral limit is 
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(a) 



(b) 



Figure 6: Two examples of pion loop graphs contributing to inr scattering in 
the quark-flow diagram picture (all lines are quark lines). Diagram (a) does not 
contain quark loops, whereas diagram (b) does. 



modified. In the quenched case the one-loop contribution gives 

(r 2 ) q y = -§| , (5.18) 

7T 

which stays finite in the chiral limit. The situation changes at two loops and 
higher: the graphs that we have discussed above, which have chiral dimension 
zero with respect to the tree level (like the two-loop fish diagram), do generate 
power-like divergences in the chiral limit. At two loops we are going to have a 
behaviour like: 

w -u- ~ j^kw (* (J W + *$ to M -) • (5 ' 19) 

where presumably also at this order the chiral logs could be reabsorbed in the 
renormalization of some order p 4 constants. 



5.4 The 7T7T scattering amplitude 

The 7T7T scattering amplitude is another example of an observable where one can 
find all the typical effects of quenching. Moreover it is an interesting quantity by 
itself since a comparison of the prediction for the two 5-wave scattering lengths 
with existing lattice calculations |l7j is possible. 



The presence of "standard" chiral logs even in the quenched theory has to be 
interpreted as due to diagrams with pion loops that do not contain quark loops. 
For the ttit scattering amplitude an example is given in Fig. ^. The one-loop 
contributions to the tttt scattering amplitude in quenched CHPT come from the 
diagrams shown in Fig. [7]. The scattering amplitude at tree level is the same as 
in standard CHPT 

A*~(s,t,u) = ^ r - , (5.20) 

where M and F are the bare pion mass and decay constant. The renormalized 
scattering amplitude in quenched CHPT and in the two degenerate flavour case 
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t-,u-channel 






s-channel 



s-channel 



Vl,2 




s-channel 



Figure 7: One-loop diagrams in quenched CHPT which contribute to the tttt 
scattering amplitude in the two degenerate flavour case. They are the "stan- 
dard" meson loop diagrams (first line) to which one has to add the corresponding 
fermionic ghost loop diagrams, the singlet vertex (x) insertion diagrams (second 
line) and the diagrams with one vi, i>2 vertex insertion (third line). 

can be written as follows 



A(s, t, u) 




+ B(s,t,u) + C(s,t,u) + 0(p 6 ) , 



(5.21) 



where 



B(s, t, u) 




aMl) 2 
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C{s,t,u) = T ^^| S 2 (2/? + /1-3) + (t-n) 2 (/l-l) 
+8sM*[l-ll + !l(Jl-l) 

+8M 4 W [JI-1 + 7 1 (If - 1) - 2 7 J (11 - l)j | . (5.22) 

For a definition of the functions J(q 2 ), h(q 2 ) and h(q 2 ) see Appendix y. The 
functions and ^(s) arise from diagrams with one and two Wq insertions on 
the two internal singlet lines in the s-channel respectively (see Fig. 0). Note that 
everything is expressed in terms of the renormalized squared pion mass M 2 given 
by Eq. ( |5.5| ) and F n = F. Note also that any dependence upon quenched chiral 
logarithms has been again reabsorbed in the B parameter contained in the renor- 
malized pion mass, as expected. The function C(s, t, u) contains only polynomial 
contributions, while the invariant function B(s,t,u) is the quenched analogue 
of the unitarity correction to the scattering amplitude in ordinary CHPT. It is 
important to note that unitarity is destroyed by the quenched approximation: 
the structure of the cuts in the one-loop amplitude is not related via unitarity to 
the real part of the tree level amplitude. Moreover, one can easily verify that the 
Fermi- Watson theorem, which relates, e.g., the imaginary part of the vector and 
scalar form factors to those of the corresponding partial waves of tttt scattering, 
is not valid in this case. 

In this particular example the violation of unitarity is also immediately seen 
in the presence of the finite functions I\{q 2 ) and ^{q 2 )-, which are not generated 
in ordinary CHPT. They have a nonzero imaginary part for s > 4M 2 that has a 
singularity at s = 4M 2 (of the type (s — 4M 2 )~ 1 / 2 and (s — 4M 2 )" 3//2 , respectively, 
see Appendix 0), which is a pure quenching artifact. These singularities have 
been already identified in [|1^, [18|]. Here we have rederived them in the a / 
case and inserted in the complete formula for the amplitude. 

Interesting quantities to be extracted from the hit scattering amplitude are 
the S'-wave scattering lengths. In Ref. [|K|] we calculated the coefficients of the 



chiral logarithms which arise in the quenched case and made the comparison with 
standard CHPT. Here we give the complete expression of the S'-wave scattering 
lengths in the isospin I = 0, 2 channels to one loop and comment on the anoma- 
lous behaviour of the isospin amplitude in the 1 = channel (which was already 
remarked in Ref. JT2|). The / = 0, 2 amplitudes are expressed in terms of the 
invariant amplitude A(s,t,u) as follows 

T°(s,t) = 3A(s,t,u) + A(t,u,s) + A(u,s,t) , 

T 2 {s,t) = A(t,u,s)+A(u,s,t) . (5.23) 

The pion scattering lengths a\ for a given isospin / and angular momentum I are 
defined by the behaviour of the partial wave amplitudes near threshold 

Ret I l (s)=q 2l {a I l +q 2 b l I + 0(q*)} , (5.24) 
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which enter the expansion in partial waves of the isospin amplitude 



r J (s,i) = 32tt^(2/+ l)P,(cos 0)tf(s) . (5.25) 
1=0 

For more details about the notation we refer the reader to Ref. || . The scattering 
amplitude in the 1 = channel contains the amplitude A(s,t,u) and therefore 
in the quenched case acquires a sick threshold behaviour due to the presence 
of functions Ii(s) and h(s). These functions do not contribute to the 1 = 2 
amplitude. On the other hand the divergences at threshold present in the infinite 
volume case show up as "enhanced" finite volume corrections to the Luscher 



formula fll9fl , that is used on the lattice to extract the scattering lengths; these 



finite volume corrections have been studied in Ref. | l2"f . We can formally define 
the quenched 1 = S-wave scattering length a[j as the coefficient of the (q 2 ) 
term in the expansion of the real part of the isospin amplitude T°(s, t) in partial 
waves. This gives us an idea of the size of normal one-loop corrections to the 
scattering length. The present definition is also equivalent to the one adopted in 



Ref. [E| in the analysis of the finite volume corrections. The quenched S-wave 
scattering length in the 1 = 2 channel a 2 , is defined in the usual way. For the 
complete renormalized S'-wave "quenched scattering lengths" at one loop we find: 

32?rF 2 _ Ml 



Ml a ° " ' 16tt 2 F2 



7 + ttt^ {7 + 5(1? + 211) + ll (511 + 1) + 2 T f (Jl - 1) - 48 Vl } 



(ml - aMl) (2 \ 5 (m 2 ~ f , 9fi x 



48vr 2 F 2 V3 J 9 A&^MlFl 



32ttF 2 , _ M\ 



-a 2 



Ml ^ 167T 2 F 2 



2 + ( 2([? + 2 ^ 2 - 1} + 2 ^ ^ ~ 1} - 47 " - 1} ) 



(m 2 - aMl) (2_ a _\ + 2 (m 2 - aMlf ^ 



24vr 2 F 2 V3 / 9 48tt 2 M 2 F 2 

The renormalized quenched scattering lengths depend upon four counterterms 
If, . . .11 and the parameters of the anomalous singlet sector at leading order, ml, 
a, vi and v-i. The counterterms l\ carry the chiral logarithms If — — logm + . . .. 
In ordinary CHPT the chiral logarithms are largely dominant in the one-loop 
corrections to the S'-wave scattering lengths at the renormalization scale /x = 1 
GeV [[2(J . Here the main unknown is the value of the parameters v± , v<i of the 



singlet sector. The singlet parameters m and a can be extracted from lattice 



calculations. Favoured values are listed e.g. in Ref. [21]| . With these values 
at hand we can do the following numerical exercise. Let us disregard for the 
moment the parameters V\ and v% and limit the analysis to the contributions 
that are reasonably expected to be the dominant ones: 1) the singlet corrections 
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in mo and a and 2) the standard chiral logarithms. With the definitions 

^ - M ' ?=j- a 5=f-«, (5.2 f 



48tt 2 F 2 ' 48tt 2 F 2 ' 48vr 2 F 2 

the leading contributions to the scattering lengths are as follows: 

7-(|«-i)» + £-'™. f + 

32ttF 2 , /2 \ = 2<5 2 , M 2 



a » = - 2 + U Q - 1 J M+ 97- 12elD / + '" • (5 - 29) 

For the numerical calculations we use F„ = 93 MeV, 5 = 0.15 and a = 0.6 
and vary the pion mass between its physical value M n = 140 MeV, and = 
600 MeV, which is presumably already outside a reasonable range of validity 
for ordinary CHPT. The chiral log is evaluated at fi — 1 GeV. The numerical 
results are given in Tables |2] and |3] for the 1 = and 1 = 2 scattering lengths 
respectively. We note that at the physical value of the pion mass the S 2 /e term 
is largely dominant in both cases: the divergence in the chiral limit produced 
by quenching is already felt at the physical pion mass. This also means that 
the whole framework is not very reliable in this range, since also higher loop 
effects may produce modifications of the same chiral order (higher powers of 5 
with the same 1/e in front). At larger values of the pion mass, which are those 
typically used in lattice calculations, the situation changes and the standard chiral 
logarithms become dominant, as it happens in standard CHPT. 













M w (MeV) 


tree 


5 




elnM 2 


140 


7 


0.09 


2.5 


1.2 


300 


7 


0.08 


0.47 


3.5 


600 


7 


0.06 


0.06 


5.9 



Table 2: Numerical values of the leading contributions to afj quenched up to one 
loop for M w = 140, 300, 600 MeV, according to Eq. flCTp . 



This picture, although at a semiquantitative level, suggests that quenched 
lattice calculations of the S-wave scattering lengths with a moderately high pion 
mass (like the ones in Ref. |I7|), should not be too far from those predicted by 
full CHPT. This conclusion is based on two observations: first the standard chiral 
logarithms start soon to be dominant with respect to the dangerous quenching 
effects, and second their coefficient happens not to be substantially changed by 
quenching ||10|| . The comparison between the standard CHPT prediction at one 
[[| and two loops |2"2" |, and the lattice calculation \T7\, has been made in Ref. |]23|| . 
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M w (MeV) 


tree 


5 


5 2 /e 


elnM% 


140 


-2 


-0.18 


1.0 


0.23 


300 


-2 


-0.16 


0.19 


0.63 


600 


-2 


-0.12 


0.02 


1.1 



Table 3: Numerical values of the leading contributions to quenched up to one 
loop for M n = 140, 300, 600 MeV according to Eq. flCTQ 

6 Summary and conclusions 

In this paper we have analyzed the quenched version of Chiral Perturbation The- 
ory at the one-loop level. We have calculated the one-loop ultraviolet divergences 
of the theory at the level of the generating functional, and shown how one can re- 
absorb all those divergences by a proper definition of the counterterms. We have 
shown that even in the presence of the anomalous singlet sector the ultraviolet 
divergent part of the quenched generating functional can be calculated in closed 
form. We have closely followed the notation and methods of standard CHPT || 
in order to identify as clearly as possible the changes produced by the quenched 
approximation in the formulation of the effective theory. 



We have found a systematic cancellation of the flavour-number dependent 
terms inside the divergent part of the generating functional to one loop. As we 
anticipated in Ref. [10 the complete iV-independence of quenched CHPT is wel- 
come, since it shows that we understand the differences between standard CHPT 
and its quenched version. Let us recall that the calculation of the divergences to 
one loop in CHPT produces explicit N dependence, in three different powers: N, 
1/N and 1/N 2 . The terms linear in N must be generated at the quark level by 
virtual quark loops: therefore they must be absent in the quenched theory. The 
terms with inverse powers of iV are generated by the decoupling of the singlet 
field from the octet of the Goldstone bosons. Since the decoupling does not take 
place in the quenched theory, also the inverse powers of N disappear in qCHPT 
to one loop. A posteriori one could say that the changes that lead from standard 
CHPT to its quenched version could have been guessed by simply looking at the 
N dependence of the generating functional to one loop. In fact this can still be 
done in other sectors of the effective theory that have not been fully analyzed 
yet. We give one example of this in Appendix [Bj, where we study the one-loop 
divergences in the sector of the on-shell non-leptonic weak interactions. 
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The quenched approximation produces a double pole in the singlet two-point 
function, which however is not allowed in a consistent quantum field theory, and is 
therefore the source of many sicknesses of quenched CHPT. As was shown already 
in Refs. [@, one of the consequences of this double pole is the appearance 
of a new kind of chiral logarithms in the one-loop corrections. Together with 
the standard M 2 In M 2 chiral logarithms, qCHPT has corrections of the form 
m^lnM 2 , which diverge in the chiral limit. The complete calculation of all the 
ultraviolet one-loop divergences in the generating functional, has shown that the 
quenched chiral logs can be accounted for via a renormalization of the lowest order 
constant B (which is proportional to the quark condensate). As a consequence, 
the renormalized B parameter diverges in the chiral limit, while the renormalized 
pion mass M 2 = 2B m q does not. The use of the renormalized pion mass to 
express any other observable makes the quenched chiral logs disappear at one 
loop, with the only exception of qq matrix elements, that are proportional to the 
renormalized B parameter alone. Hence, qq matrix elements remain the unique 
place for discovering the presence of quenched chiral logs in quenched lattice 
calculations within the strong sector. 



The double pole in the singlet two-point function also changes the standard 
chiral power counting for which diagrams with a higher number of loops are of 
higher chiral order. In the quenched case one may have graphs with any number 
of loops with the same chiral dimension, and the chiral order of an amplitude 
is no more constrained to be positive: as a consequence, quenched CHPT has 
power-like divergences in the chiral limit. These divergences are in principle 
a very serious problem of the theory, although they seem to be a unavoidable 
consequence of the quenched approximation. Since the graphs that have negative 
chiral dimension are also ultraviolet finite, their study requires the calculation of 
the UV finite part of the loop corrections. We have shown how they arise within 
the generating functional approach; at one-loop and at order $ 4 they are given 
in Eq. ( |3.22| ). We have therefore analyzed some physical quantities at one loop 
in the case of two degenerate light flavours: the scalar quark condensate, the 
pion mass, the scalar and vector form factors of the pion and the scattering 
amplitude. This has given us the possibility to discuss in detail the changes 
induced by quenching in the UV finite part of the one-loop corrections. The main 
changes can be summarized by saying that unitarity is not satisfied anymore, and 
that the double pole in the singlet two-point function produces singularities in 
the chiral limit, and also unphysical singularities at threshold. 



The differences between CHPT and its quenched version are rather well under- 
stood, as the study of the flavour-number dependence of the generating functional 
at one-loop also shows. The presence of the double pole in the singlet two-point 
function is also a rather direct consequence of the quenched approximation. This 
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double pole has dramatic effects on the effective theory. However, it looks plau- 
sible that despite all these inconsistencies (or maybe because of them) quenched 
CHPT is the right tool to understand the effects of quenching in actual lattice 
calculations. The crucial check will be a detailed comparison of qCHPT pre- 
dictions with the quark mass dependence of various quenched quantities on the 
lattice, and especially of the way they approach the chiral limit. We expect that 
further investigations in this direction will answer these questions. 
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A Divergences proportional to m§ 

In this appendix we try to explain in a simple way why the divergence pro- 
portional to wig is given by the chiral invariant term (x+)- We start from the 
observation that the double pole in the singlet propagator can be expressed with 
a derivative of a normal propagator: 

ml m 2 9 1 (A 1) 



(M 2 -p 2 ) 2 dM 2 M 2 -p 2 



In fact in the quenched case one is keeping only the first two terms in a Taylor 
series: 

(MZ-M 2 ) —— - +... , (A.2) 



M 2 -p 2 M 2 -p 2 v * ' (M 2 - p 2 

where m\ = M 2 , —M 2 . The divergences proportional to m\ arise from this Taylor 
expansion of the propagator. Before the expansion, singlet loops give two type 
of divergent contributions: 

1 f 1 1 f d d p 1 1 , K rj\ 

and 7 / 77777772 77777 77 77 ■ ( A - 3 ) 



% J (2n) d M 2 - p 2 ' i J (2n) d M 2 - p 2 M 2 - {p - q) 

The divergence proportional to m\ is obtained expanding the singlet propagator 
inside the loop integrals and taking only the second term in the expansion. It is 
clear that only the tadpole produces a divergence: 

I r d d p 1 1 r d d p -m 2 _ 2 

iJ [2-kYM 2 -p 2 i J (2nY(M 2 -p 2 ) 2 ~ [ J ' 1 J 

The chiral structure of the term proportional to the tadpole is very easily iden- 
tified, and is given by the term proportional to £q after expanding the action in 
fluctuations around the classical solution: 

C 2 = ^(uX + X+ )-^{-^(o + M 2 )^-e ((x + )-2NM 2 )} . (A.5) 

Obviously, the tadpole is generated by the contraction of the two £ ' s in the last 
term of Eq. (|A.5[) , and therefore the divergence is proportional to (x+)- 

In summary: the divergence proportional to ml must be, for dimensional 
reasons, a chiral term of order p 2 . This divergence comes out from a tadpole 
through a derivative with respect to M 2 ,. In the tadpole the vertex is a chiral 
invariant of order p 2 . The simplified recipe to derive the form of the divergent 
term proportional to ml amounts to determining the chiral invariant that is 
proportional to £g after expanding the CHPT Lagrangian in fluctuations around 
the classical solution. 

The recipe applies in the same way to other sectors. We will give an example 
below for the non-leptonic weak interactions sector. 
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B Non— leptonic weak interactions 



In this appendix we show how one can derive the structure of the divergences 
in the quenched case just by looking at the N dependence of the divergences 
of the standard CHPT case. We will consider the weak octet Lagrangian which 
contributes to the non-leptonic weak interactions with AI = 1/2. The structure 
of the divergences in standard CHPT has been given by Kambor, Missimer and 
Wyler and then expressed in terms of a minimal basis by Ecker, Kambor 
and Wyler [^5[] for on-shell processes. We will use the basis given in the latter 
Reference. 

Let us recall here some basic notation. In this sector the lowest order La- 
grangian is given by: 

4 2 = c 2 (A«/)+c 5 (Ax + ) , (B.l) 

where A = u\u\ and 02,5 are low energy constants. The C5 term can be omitted 
for on-shell processes, as it can be transformed away by a field redefinition. We 
will not consider it anymore in what follows. Since the above-mentioned calcu- 
lations of the divergences were made for N = 3 we have redone the calculation of 
the divergences for N generic. Our result for the divergent part of the one-loop 
generating functional, reads as follows: 

,d-4 



IT 



C>2 



one loop 



16tt 2 d - 4 Fl 



dx L A . + finite terms 

div 1 



W4 + \ws - ~W 7 + ^W 8 - ±W n 



+ 



1 

N 



+ N 



2W l0 + Wi 2 - 2W 21 - 2W 22 + W 36 - —Wu 
\w x - \w 2 + \w h + ]w 9 - V 12 + -Wu + V15 - ^w 16 

3 o 4 4 4 12 o 12 



1 



— w w 

12 9 



Wis 

24 8 

1 1 1 

_i — w — iy 36 — w 37 

24 4 24 



-W20 + \~W21 + \w 22 + l\V 25 
4 2 2 o 



4 



-W, 



:-!S 



(B.2) 



With Wi we have indicated the operators of order p A given in Ref . ] . Here we 
have one more: 

W 38 = (Au^x+u* 1 ) , (B.3) 

that for = 3 is linearly dependent on the other operators thanks to the trace 
identity coming from the Cayley-Hamilton theorem: = —W5 + W6 + I/2 W 7 + 

According to the rules we have found in the strong sector, the divergences of 
the octet weak sector for on-shell processes in the quenched case become: 



div 



(B.4) 
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For the sake of clarity we list the definitions of the W{ terms needed here: 



W, = (AuJ {u»u v u v ) , W 6 = (A Ufi ) (x+u^) , W 7 = (A X+ ) > 
J^ 8 = (A«X) , = (Ax+) (x+> • (B.5) 

We have verified that the contribution of the singlet sector to the divergences in- 
deed cancels the terms with a negative power of N. It would remain to be checked 
by an explicit calculation that the inclusion of the fermionic ghosts removes the 
terms with the linear flavour dependence |L3||, as one expects. 



Finally, we derive the contribution of quenched chiral logs in this sector, ac- 
cording to the simple recipe given in the previous appendix. The answer is very 
simple: there are no quenched chiral logarithms in the on-shell octet weak sector 
at one loop. The reason is that if one expands the term proportional to C2 in the 
lowest order Lagrangian around the background field, there are no terms with 
£q. The situation would be different considering also off-shell processes, since the 
expansion of the C5 weak mass term generates the singlet term cs£o(A% + ) . The 
conclusion here is that the quenched chiral logs can be reabsorbed in a redefinition 
of c 5 . 



C One— loop functions 

Here we list the functions J(q 2 ), Ii(q 2 ), ^{q 2 ) which appear in the one-loop cor- 
rections to the quenched observables analyzed in this paper. h(q 2 ) and hiq 2 ) 
are not generated in standard CHPT and arise from the insertion of the a and 
m 2 , vertices in any internal singlet line. 

The one-loop function J(q 2 ) in Minkowski space-time, is given by 

^ q ' = ~il (2vr) d (M 2 - P) (M 2 -(I- qf) ' ^ ' 

and 

J(q 2 ) = J(0) + J(q 2 ) , (C.2) 
where J(0) contains the divergent part 

1 / M 2 \ 
J(0) = -2A-— - In — + 1 +0(d-4) , 



A 



/i 



16tt 2 

d-4 r 1 1 



16tt 2 



d-A 2 



-(ln47r + r'(l) + l) 



(C.3) 



and J{q 2 ) is finite. The explicit expression of J(q 2 ) for d — 4 is: 



a 



J ^ = T^-{^—i + 2 > • (a4) 
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where a = sfl - 4M 2 /q 2 . The UV finite functions h{q 2 ) and J 2 (g 2 ) are given by 

d d l 1 



W) = ■ 



(2vr) d (M 2 - Z 2 ) 2 (M 2 - (Z - g) 2 ) 
d d l 1 



(2vr) d (M 2 - I 2 ) 2 (M 2 - (/ - g) 2 ) 2 
They are the F.T. of the functions h{z) and hiz): 

h(q 2 ) = -Jd d ze^ z h{z) , 

h(q 2 ) = - I d d ze iqz I 2 {z) . 
% J 

The explicit expression of h(q 2 ) and hiq 2 ) for d — 4 is as follows 



1 11 (7-1 
167T 2 g 2 (T (7+1 

1 1 



, g 2 - 2M 2 , a - 1 
1 + In 

(7+1 



g 2 a 



8tt 2 gV 2 

In the text we have used also these functions subtracted at g 2 = 0: 

/i,2(g 2 ) = /i,2(0) + / li2 (g 2 ) , 

with 



h(0) 



1 



,/ 2 (o) = 



1 



(C.5) 



(C.6) 



(C.7) 



(C.8) 



32vr 2 M 2 '^ v " y 96tt 2 M 4 ' 
When g 2 > 4M 2 both functions develop an imaginary part, which is given by 

1 1 1 



(C.9) 



lmh(q 2 ) = 
ImJ 2 (g 2 ) = 



167T g 2 a ' 

1 2(g 2 -2M 2 ) 



16tt 



g 3 <7 3 



(CIO) 



Notice that they diverge at g 2 = 4M 2 
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